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ABSTRACT 

We present a proof-of-concept study that dust extinction curves can be extracted from the infrared 
(IR), optical, ultraviolet (UV) and X-ray afterglow observations of GRBs without assuming known 
extinction laws. We focus on GRB 050525A (z — 0.606), for which we also present IR observations 
from the Spitzer Space Telescope at t = £ir w 2.3 days post-burst. We construct the spectral energy 
distribution (SED) of the afterglow at t = iiR and use it to derive the dust extinction curve of the host 
galaxy in 7 optical/UV wavebands. By comparing our derived extinction curve to known templates, 
we see that the Galactic or Milky Way extinction laws are disfavored versus the Small and Large 
Magellanic Cloud (SMC and LMC) ones, but that we cannot rule out the presence of a LMC-like 
2175 A bump in our extinction curve. The dust-to-gas ratio present within the host galaxy of GRB 
050525A is similar to that found in the LMC, while about 10 to 40% more dust is required if the 
SMC template is assumed. Our method is useful to observatories that are capable of simultaneously 
observing GRB afterglows in multiple wavebands from the IR to the X-ray. 

Subject headings: gamma rays: bursts — galaxies: ISM — dust 



1. INTRODUCTION 

Extinction 8 of light plays a major role in the obser- 
vation of any astronomical object, but even more so for 
intermediate- and high-redshift ones due to the intrin- 
sically fainter signals and our poorer understanding of 
the dust physics involved. In the optical range of wave- 
lengths, dust grains are the dominant contributor to ex- 
tinction. High-quality extinction curves have been con- 
structed for our Galaxy (the Milky Way; MW), as well 
as for the Small and Large Magellanic Clouds (SMC and 
LMC; e.g., Cardelli, Clayton & Mathis 1989; Pei 1992; 
Weingartner & Draine 2001). 

Despite the importance of the subject, surprisingly lit- 
tle is presently known about the extinction properties 
of dust grains in the universe at redshifts of z > 0. 
Some constraints on the properties of dust have been 
obtained for active galactic nuclei (AGN; Laor & Draine 
1993) and starburst galaxies at low redshifts (Calzetti 
1997). Nevertheless, our knowledge of the composition 
and grain size distribution of dust in intermediate- and 
high-z galaxies is very limited. The only constraints set 
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on some intermediate- z galaxies (z ~ 0.5 to 0.9) were 
obtained through the analysis of multi-imaged quasars 
(QSOs) behind a lensing galaxy (Falco et al. 1999; Motta 
et al. 2002; Munoz et al. 2004). These studies utilize the 
fact that the lensed images have the same intrinsic spec- 
trum. The extinction curves, which are measured and 
characterized by a Cardelli, Clayton & Mathis (1989) 
profile, range from flat ("gray") to steep, SMC-like pro- 
files, both with and without the 2175A feature. 

Other studies have derived attenuation curves, ob- 
tained by substituting a source population — for which 
an average, intrinsic spectrum can be safely assumed — 
for the point-like source. Notable results were achieved 
by Calzetti (1997) for starburst galaxies and by Johnson 
et al. (2007) for a sample of 1000 Sloan Digital Sky Sur- 
vey (SDSS) galaxies. Strikingly, the former and latter 
studies yielded gray and SMC-like attenuation (or "ef- 
fective extinction") curves, respectively. These studies 
leave open the question of whether the inferred diversity 
stems from intrinsic differences in the dust properties or 
is an effect of hidden (and incorrect) assumptions in the 
utilized methods. 

Our relatively poor knowledge of dust physics at z > 
makes it difficult to derive the fundamental properties of 
intermediate- and high-z galaxies such as the star forma- 
tion rate (or at least its obscured fraction), the metallic- 
ity of the interstellar medium (ISM), the rate of super- 
nova (SN) explosions, and the extinction-corrected lumi- 
nosity of any galactic source of light. Our level of under- 
standing of dust physics in our Galaxy shows us what can 
be learned if high-quality extinction curves are available. 
Studying the extinction properties of the ISM in the MW 
has provided insight into dust composition, which has in 
turn taught us about grain coagulation and destruction, 
metal depletion (into dust), and grain polarization prop- 
erties (e.g., Mathis, Rumpl & Nordsieck 1977; Weingart- 
ner & Draine 2001; Draine & Li 2007). These pioneering 
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studies have set the stage for asking whether MW dust 
physics is unique to our Galaxy or evolves with redshift. 
Differences between the MW, SMC and LMC extinction 
laws suggest that grain properties probably vary with 
the local radiation field, metallicity and galactic mass, 
but the lack of more examples has hampered any further 
conclusions. 

Measuring the extinction curve at intermediate- or 
high-z is a formidable task, due to the rarity of sources 
with known spectra and luminosities. The discovery 
of 7-ray bursts (GRBs) and their afterglows, about a 
decade ago, has opened a new window into the study 
of dust physics. In particular, long-duration GRBs are 
believed to occur in dusty environments, since it is 
thought that they result from the deaths of very mas- 
sive stars. Their afterglows have featureless, power-law 
spectra (e.g., Liang et al. 2007) and intense luminosi- 
ties (~ 10 48 to 10 50 erg s _1 ), easily outshining a typical 
QSO (~ 10 46 erg s" 1 ). It has been previously demon- 
strated that the study of the infrared (IR), optical, ul- 
traviolet (UV) and X-ray afterglows of a GRB can place 
constraints on the extinction by dust present within the 
host galaxy (e.g., Price et al. 2002; Starling et al. 2007). 
However, these studies lacked an unbiased measurement 
of the unabsorbed spectrum and its normalization, re- 
stricting the analyses to only fitting known extinction 
laws to the data, rather than directly measure the dust 
extinction. 

In this paper, we pioneer a method that is similar in 
spirit to the previous studies, but uses less assumptions. 
For our study, we focus on GRB 050525A, the first detec- 
tion of a GRB at mid-IR wavelengths with the Spitzer 
Space Telescope (Garnavich et al. 2005). We demon- 
strate that by constructing a spectral energy distribution 
(SED) of the GRB afterglow which includes IR and X- 
ray data points, one can derive an extinction curve with- 
out having to assume a known dust extinction law. In 
other words, our method allows a template-independent 
calibration of the extinction curve. We emphasize that 
such an approach is well-suited for observatories that are 
capable of simultaneously observing GRB afterglows in 
multiple wavebands (e.g., GROND). 

In SJ21 we describe the IR, optical, UV and X-ray obser- 
vations used in our study. Our Spitzer IR observations at 
t = <ir « 2.3 days are not contemporaneous with the op- 
tical, UV and X-ray observations used, which motivates 
a description, in Sj3l of our light curve fitting procedure 
that interpolates for the optical, UV and X-ray SED data 
points at t = Using our constructed SED, we derive 
an extinction curve for the host galaxy of GRB 050525A 
in §3J we further analyze its dust-to-gas ratio if known 
extinction laws are fitted to our curve. Finally we com- 
pare our results to that of previous studies, discuss their 
implications and present opportunities for future work in 

m 

2. OBSERVATIONS 

2.1. Optical, Ultraviolet & X-ray Data from Literature 

GRB 050525A occurred at 00:02:53 UT on 2005 May 
25th (Band et al. 2005) at a redshift of z = 0.606 (Fo- 
ley et al. 2005). Optical and UV observations of its 
afterglow, in the UVW2 (198 nm), UVM 2 (220 nm), 
UVWl (260 nm), U, B and V bands were taken by 
Blustin et al. (2006, hereafter B06) using the Swift 



UVOT. B06 also performed Swift XRT measurements 
in the 2 to 10 keV bands. The optical and UV obser- 
vations were corrected for Galactic extinction using the 
extinction data of Schlegel, Finkbeiner & Davis (1998) 
and Pei (1992). For corrections to the X-ray observa- 
tions, the Galactic hydrogen column density was taken 

to be 2Vg = 9.0 x 10 20 cm~ 2 ; the hydrogen-equivalent 
column density at z — 0.606 was estimated to be about 
2—3 x 10 21 cm" 2 . (See §2.3 of B06.) 

The R band data is taken predominantly from Klotz 
et al. (2005) and Delia Valle et al. (2006, hereafter D06) 
for early (< 100 min) and late (> 100 min) times, re- 
spectively. These are supplemented by the observations 
of Cobb & Bailyn (2005), Homewood et al. (2005) and 
Haislip et al. (2005). GRB 050525A is thought to be 
observed on-axis and is associated with supernova (SN) 
2005nc; D06 proposed that a bump in the R band light 
curve at ~ 20 days is due to a SN contribution peak- 
ing at between 10 and 20 days. We therefore consider R 
band data only up to about 5 days. The R band magni- 
tudes are converted to flux densities taking into account 

the Galactic extinction of AjP = 0.25 mag towards GRB 
050525 A (Schlegel, Finkbeiner & Davis 1998). 

To account for systematic errors inherent in using data 
from different observatories, we enforce a minimum error 
of 5% (~ 0.05 mag) for the optical/UV data. 

2.2. Infrared Observations from Spitzer 

Our IR observations were triggered using a "target-of- 
opportunity" (ToO) program on the Spitzer Space Tele- 
scope (Garnavich et al. 2005). Imaging and spectroscopy 
began on May 27.264 (UT) at t = t m = 1.95 x 10 5 s w 2.3 
days post-burst. We obtained measurements (Table [l| 
from the 3.6, 4.5, 5.8 and 8.0 fim IRAC channels, and 
from the 24.0 /jm MIPS channel, using the latest Spitzer 
Science Center pipeline products (S16; 2007 Dec). The 
most current flux calibration was used (Engelbracht et 
al. 2007). The data required very little post-processing 
beyond the standard products, because the source was 
placed at the center of the array (Engelbracht et al. 
2007). An aperture photometry measurement was car- 
ried out using a 7"radius at the position of the source 
(RA = 18h32m32.58s, DEC = 26d20m22.5s, J2000) that 
required an aperture correction of 2.05 (Table 1 in En- 
gelbracht et al. 2007) to recover all of the flux. The 
selected aperture radius reduces the light contamination 
by a nearby source; for the same reason, the background 
was measured on a comparable area ~ 10" away from 
GRB 050525A. 

2.3. Host Galaxy Contribution & Extinction 

The optical and ultraviolet contributions from the host 
galaxy are at the < 1% level on average, especially at 
early times (B06). There is a definitive measurement 
of the host galaxy contribution only in the R band: 
R = 25.2 mag (D06). Since the X-ray and other opti- 
cal/UV contributions are expected to be small, we ignore 
them in our analysis. However, it is conceivable that the 
host galaxy and GRB afterglow emissions might be of 
comparable strengths in the 24.0 /im band (e.g., Zheng 
et al. 2006; Martin et al. 2007). Unfortunately, the host 
galaxy contribution at 24.0 /xm cannot be empirically 
determined with just a single epoch of IR observations. 
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We therefore consider two bracketing scenarios: Models 
A and B. In Model A, the measured 24.0 /xm flux is as- 
sumed to be entirely from the GRB afterglow. In Model 
B, we reduce the 24.0 /Ltm flux measurement by a fac- 
tor of 2 when constructing a spectral energy distribution 
(SED) at t = tm « 2.3 days post-burst (see g}. 

We further assume that dust extinction by the host 
galaxy in the IR and X-ray is negligible. We show later 
(in M4. 3[) that the derived visual extinction is Ay ~ 0.2, 
which implies, for example, that dust extinction in the 
3.6 /im band is negligibly small, since it is ~ 5% of that 
in the V band (B. Schmidt 2008, private communica- 
tion). However, there may be a gray component to the 
attenuation (e.g., Perley et al. 2008), which we do not 
consider. 

3. LIGHT CURVE FITTING 

Due to the lack of optical and UV data taken at the 
same time as the IR observations, we evaluate the fluxes 
at t = tm rs 2.3 days by modeling the light curves ac- 
cording to the standard afterglow model (Sari, Piran & 
Narayan 1998; Panaitescu & Kumar 2000). The thor- 
ough sampling of the R band light curve at late times 
allows for an interpolation, and not a much more dan- 
gerous extrapolation, of the data. We assume that the 
optical and UV radiation is powered by the same popula- 
tion of electrons, thus allowing for the temporal evolution 
of the optical and UV light curves to be the same. 

We fit to the data a double power law of the form 

tp li, A - TP li, i \ [" (W^)^ 1 + (tm./t b ) sa ~ : 
v{ ' >- vK , m) [ (t/t b ) sai + {t/t b ) sa2 \ ' 

where ot\ and oti are the asymptotic power law indices 
at t <C tb and t ^> %, respectively, and v is the observed 
frequency. The index s controls the sharpness of the 
transition at t = ty, we adopt s — 3 throughout, since 
it turns out to be irrelevant for our results (see Liang et 
al. 2007). The transition at t = i& may be interpreted 
as a jet break; its existence at ~ 0.16 to 0.2 days has 
been discussed in the literature (B06; Liang et al. 2007), 
albeit with some controversy (Sato et al. 2007). D06 
reported a jet break at about 0.3 days from an analysis 
of the R band data only. 

For the optical and UV data, we fit for F u {u, tm), ai, 
a.2 and tj, in the 7 wavebands simultaneously (i.e., 10 
parameters; Fig. [JJ). The optical/UV fit is largely con- 
strained by the R band, for which we have the most data 
points (44). The reduced chi-square value for the com- 
bined optical/UV fit is xl = 310/142 = 2.2. Since the 
excess \ 2 is evenly distributed at all times and frequen- 
cies, we believe this somewhat large value of x r to be 
due to either an under-estimation of the data errors or 
to intrinsic high-frequency fluctuations in the light curve. 
This will be appropriately taken into account when com- 
puting the uncertainties on the fitted parameters. 

We perform an analogous fit to the X-ray data, which 
yields x T = 68.5/54 = 1.3. Upon obtaining a model fit 

(x) 

to the X-ray flux, we convert it into a flux density, P; . 
Following B06, we assume a power law with a spectral 
index (3 X w -0.9, Pj x) oc iA. 

The evaluation of the uncertainties in the optical and 
UV fluxes at t = tm is of particular importance for our 



study. For each parameter P, we compute the marginal- 
ized chi-square, Xm(P)- The 1-a errors, P CT , are deter- 
mined from the condition: 

X 2 m (P.) - mm{ X 2 m (P)} = Xr- (2) 

The preceding condition is stronger than the usual one 
(where the difference between the LHS quantities is set 
to unity), because it takes into account the uncertainties 
due to random fluctuations in the data that cause the 
minimum % 2 to have a higher than optimal value. The 
resulting fluxes and their uncertainties at t = £tr are 
reported in Table [TJ 

For the optical and UV data, we obtain: ai )OV = 0.95± 
0.01, a 2 , ov = 1.80±g;gg and i h , ov = (2.50±g;|?) x 10 4 s 
ks 0.3 days; our tb )0 v value is consistent with the finding 
of D06. For the X-ray data, we get: a x , x = OJOlg^, 
a 2 , x = 1.56 ± 0.04 and t M = (1.46±g;|f) x 10 3 s « 0.4 
hours. 

4. RESULTS 

4.1. Constructing the SED at 2.3 Days Post-Burst 

Using the optical, UV and X-ray fits, as well as the 
IR data and upper limit, we construct the SED for the 
afterglow at t = tjR for both Models A and B (Figs. [2] 
and [3]). In the standard fireball model of Sari, Piran & 
Narayan (1998), the SED at t = tin transitions from a 
power law with an index of (1 — p) /2 to one with an index 
of — p/2 at the cooling frequency of 

„ c = 3.8 x 10 16 Hz e-%1 (^) e£», (3) 

where 6b,-3 is the fraction of the shock energy injected 
into the magnetic field in units of 10~ 3 , no is the den- 
sity of the circum-burst material in units of 1 cm~ 3 , 
Piso = Piso,53l0 53 erg is the isotropic equivalent energy, 
andg c = (p— 0.46) exp (— l.ldp) is a function of the power 
law index p m 2.5 of the electron distribution. (Sec also 
Shao & Dai 2005.) With our SED, we do not see evi- 
dence for a cooling break, consistent with the findings of 
B06, who noted the absence of a break for the SED at 
t = 25,000 s. 

The observed absence of a synchrotron cooling break 
in the SED allows us to fit a single power law to it, us- 
ing only the IR and X-ray flux densities. The inferred 
indices are (3 = -0.87 ± 0.02 { Xr = 3.6/3 = 1.2) and 
P = -0.85±0.02 (xl = 2.6/3 = 0.9) for Models A and B, 
respectively. Uncertainties in the unabsorbed continuum 
are evaluated filter by filter; centering the frequencies at 
each filter frequency v' allows us to avoid covariances 
with the uncertainty in the SED slope. Namely, for 

In F v = f3 lnx + lnP^o, (4) 

where x = v /is' , AP^ = AP„.o when x = 1. 

4.2. Deriving the Extinction Curve 

For Models A and B, the 7 optical and UV data 
points lie below the SED fit, indicating extinction by 
dust present within the host galaxy of GRB 050525A. 
The extinction in a given waveband is related to the op- 
tical depth, t\, by A\ = 2.5T>/hilO; we calculate A\ for 
our 7 optical/UV wavebands and present them in Ta- 
ble [2j The derived dust extinction curves are shown in 
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Figs. [2] and [3] for Models A and B, respectively. In gen- 
eral, the extinction curve is a steeply-increasing function 
of the photon frequency in the rest frame of the host 
galaxy, z/o, and can be fitted with a power law: the re- 
sulting indices are 6 = 1.01 ± 0.31 {\ 2 r = 5.6/5 = 1.1) 
and S = 1.23 ±0.40 (xj = 5.5/5 = 1.1) for Models A and 
B, respectively. These fits are noticeably steeper than 
the t\ oc i/q 7 attentuation curve derived by Johnson et 
al. (2007). Broad features, such as the 2175 A bump 
seen in our Galaxy and the LMC, do not appear to be 
present (but see £|4.3|) . Starling et al. (2007) derived sim- 
ilar results for their sample of 10 GRB hosts using known 
extinction law templates. 

We note that if we only use the 3.6 and 4.5 /im flux 
measurements when fitting a power law function to the 
SED, then our results are consistent with zero extinction 
in the V and R bands. We tabulate the extinction values 
resulting from such a scenario as "Model C" in Table [H 
In the other extreme, we explore a "Model D" , where we 
instead use only the 8.0 and 24.0/^m (and ignore the other 
IR) data points, which serves as a conservative upper 
limit on the amount of extinction present . 

We caution that the Galactic extinction towards GRB 
050525A is comparable to, and in some cases larger than, 
our derived intrinsic extinction. Given the possible vari- 
ations in the Galactic extinction curve, it is possible 
that additional sources of errors are present which we 
have not considered. For example, Schlegel, Finkbeiner 
& Davis (1998) report a Galactic y-band extinction of 

= 0.316, but with an error of AA$ ~ 0.1 or ~ 30%. 
Other sources of uncertainty include the Galactic UV 
bump (Savage et al. 1985) and the assumed Ry value 
for the MW (Cardelli, Clayton & Mathis 1989; and ref- 
erences therein). 

4.3. Obtaining Dust-to-Gas Ratios for Various 
Extinction Laws 

We now wish to study the consistency of our results 
with the MW, SMC and LMC extinction laws. These 
extinction laws were previously presented in Weingart- 
ner & Drainc (2001) and Li & Draine (2001), but have 
since been updated in Draine (2003a, b, c). We first con- 
volve the model extinction laws with the relevant filter 
responses, which for simplicity are assumed to be Heav- 
iside (i.e., "top hat") functions. The filter widths are 
taken to be 128, 75, 98, 88, 70, 51 and 76 nm for the R, 
V, B, U, UVW1 (260 nm), UVM2 (220 nm) and UVW2 
(198 nm) bands, respectively. We then fit for the optimal 
value of the host galaxy hydrogen column density 9 , Ah , 
that allows the model law (A\/N-£i) to be matched to the 
observed curve {A\). We next use our power law fits to 
the extinction curves (see £|4.2[) to derive the value of the 
visual extinction in the V band, Ay , in the rest frame of 
the host galaxy: A v = 0.22 ±0.06 (Model A), 0.15 ±0.06 
(Model B), 0.02±£; 3 (Model C) and 0.41 ± 0.14 (Model 
D). The resulting Ay /Ah values for Models A and B are 
reported in Table [3J 

We first see that the MW extinction laws, regardless 
of the Ry value and the assumption of either Model A 
or B, are somewhat disfavored even allowing for different 

9 Specifically, the hydrogen-equivalent column density intrinsic 
to the host galaxy, along the line of sight to the GRB. This could 
be different if GRBs occur in special regions of galaxies. 



dust-to-gas ratios; less dust is generally required than ex- 
pected. This conclusion is also true for Models C and D. 
We also tabulate the ratio of our derived Ay /Nh value to 
the theoretically expected one for the various extinction 
laws, a quantity we denote by fa- We see that the SMC 
and LMC laws provide good fits to our extinction curve. 
Based on the derived fd values, about 10 to 40% more 
dust is present if we use a SMC fit. For the LMC fit, 
Model A requires 10% more dust, while Model B needs 
20% less dust. The similar Xr values of the SMC and 
LMC fits do not allow us to rule out the possibility that 
an LMC-like 2175 A bump may be present in our extinc- 
tion curve. We note that Models C and D favor the SMC 
law over the LMC one. 

In general, the dust-to-gas ratio found in the host 
galaxy of GRB 050525A resembles that found in the 
LMC, and is intermediate between the ratios for the SMC 
and our Galaxy. It is also similar to that found in the 
circum-nuclear region of AGNs (Maiolino et al. 2001). 
The relationship between JVh and Ay has been studied 
in GRB hosts by Galama & Wijers (2001) and Stratta 
et al. (2004). Their derived dust-to-gas ratios are ap- 
proximately ten times smaller than in the MW, which 
is roughly consistent with our findings. It is comfort- 
ing that our derived values for the extinction curve and 
Ay /N-[{ are somewhat similar for Models A and B. 

5. DISCUSSION 

5.1. Comparison to Other Studies 

Previous studies of the dust extinction laws in the 
circum-burst environment of GRBs used known extinc- 
tion templates (MW, LMC and SMC). Assuming Galac- 
tic dust conditions, Galama & Wijers (2001) found Ah ~ 
10 21 to 10 23 cm" 2 and low extinctions {Ay ~ 0.7) for 
their sample of 8 GRBs. For 13 GRB afterglows, Stratta 
et al. (2004) examined several extinction laws and in- 
ferred similar column density values that are generally 
larger than expected for the given Ay. In most cases 
steep, SMC-like extinction curves were preferred (Star- 
ling et al. 2007), although in some cases flat, gray curves 
have been claimed (e.g., Stratta et al. 2005). Savaglio, 
Fall & Fiore (2003) infer Ay « 1 mag, in contrast to 
Ay ^0.1 mag found in damped Lya (DLA) systems 
along QSO sight lines, supporting the idea that optically- 
selected QSOs mainly probe regions of high-z galaxies 
with less dust and smaller gas column densities, while 
GRBs probe the denser, star-forming, ISM regions, and 
thus the two approaches are complementary (See also 
Savaglio & Fall 2004.) 

The Spitzer IR observations of GRB 050525A allow us 
to determine the extinction curve in its host galaxy while 
making less assumptions than in previous studies. Re- 
cently, Li, Li & Wei (2007) applied a similar but rather 
risky methodology, where they attempted to determine 
the spectral slope of the unabsorbed continuum by em- 
ploying extrapolation from the X-ray to the optical range 
of wavelengths. Lacking the IR information means that 
their results are subjected to large uncertainties. In addi- 
tion, the degeneracy of the X-ray spectral slope with the 
amount of soft X-ray photo-absorption present makes the 
computation of Ay/Nn extremely challenging. In our 
treatment, the X-ray spectral information is used only 
as a consistency check. Li, Li & Wei (2007) claimed that 
the GRB hosts they study have gray extinction curves, 



Extinction Law in Intermediate- z Galaxy 



5 



but this is not surprising considering their basic assump- 
tion regarding the slopes of their SEDs. 

5.2. Implications for Host Galaxy Studies 

A detailed study of the properties of the host galaxy of 
GRB 050525A is beyond the scope of our paper, but we 
briefly discuss some implications of our results to such 
a task. The origin of the 2175 A bump is a subject 
of controversy — Draine & Li (2007) suggest that it is 
due to the presence of polycyclic aromatic hydrocarbons 
(PAHs) in the dust population. However, these molecules 
have been found to be deficient in low-metallicity, high 
radiation-density, blue dwarf galaxies, the favored hosts 
of GRBs (e.g., Hogg et al. 2005; Engelbracht et al. 
2008). If our derived extinction curve can be more closely 
sampled in frequency, it will allow us to make a more 
definitive statement regarding the absence/presence of 
the 2175 A bump, which may shed light on the PAH hy- 
pothesis. In their study of dwarf galaxies, Engelbracht et 
al. (2008) inferred dust-to-gas ratios that differ by 3 or- 
ders of magnitude, preventing a meaningful comparison 
with our derived Ay /Ah values. 

D06 analyzed the optical spectrum of GRB 050525A at 
the time of the SN bump and concluded that SN 2005nc 
is akin to SN 1998bw, an energetic Type Ic SN (Patat 
et al. 2001), but dimmed by ~ 0.3 mag. If we speculate 
that SN 2005nc is a Type Ic event as well, it is interesting 
to ask if GRB 050525A/SN 2005nc erupted in the bright- 
est region of a spiral or irregular galaxy, where conditions 
are appropriate for star formation, as suggested by Kelly, 
Kirshner & Pahre (2007), and if the inferred dust-to-gas 
content is consistent with such an environment. Such a 
determination is relevant to the hypothesis that metal- 
poor and metal-rich Wolf-Rayet stars, the proposed pro- 
genitors of Type Ic supernovae, result in long GRBs and 
less energetic core collapse explosions, respectively (Yoon 
& Langer 2005; Woosley & Heger 2006). 

5.3. Is the Derived Extinction Curve Pristine? 

The dust extinction curve at late times (i.e., during 
the afterglow phase) may be indicative of the residual 
extinction present after the onset of the burst radiation. 
The effects of dust sublimation can be dramatic if the ab- 
sorbing region is small and compact (e.g., in a molecular 
cloud) , while they are negligible for typical ISM densities 
(Waxman & Draine 2000). A direct way to quantify the 
modification of the dust grain population by the burst ra- 
diation is to measure the optical opacities at very early 
times, when these changes are occuring. However, these 
times are generally too short for the effects to be observ- 
able (Perna et al. 2003). 

If the absorbing region is small and compact, then vari- 
ations in the optical opacities should be accompanied by 
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changes in the effective column density due to the ioniza- 
tion of the gas. These effects occur on a reasonably long 
timescale and are detectable (e.g., Lazzati & Perna 2002, 
2003; Campana et al. 2007). To within the measure- 
ment uncertainties, B06 detected no changes in Ah from 
t = 250 to 25,000 s (see their Table 6). To accomodate a 
situation in which dust is destroyed without an accom- 
panying, significant ionization of the gas requires spe- 
cially fine-tuned conditions (Perna & Lazzati 2002). We 
therefore conclude that our measured extinction curve is 
pristine to the host galaxy of GRB 050525A and was not 
modified by the prompt emission of the GRB. 

5.4. Future Work & Opportunities 

Our method of deriving the dust extinction curve of a 
GRB host galaxy is comparatively robust, since we have 
cut down on the number of assumptions made in previ- 
ous studies. Nevertheless, we regard our work as a "proof 
of concept" study, which allows extinction curves to be 
derived independently of the known templates. Future 
work may improve on ours by not only obtaining more 
IR and optical/UV observations simultaneously (so as to 
avoid the fitting procedure described in and earlier 
in time, but also sampling the SED at more frequen- 
cies in between the IR and the X-ray. Such observations 
should be complemented by studies of the host galaxy 
emission. For example, the Gamma-Ray Burst Opti- 
cal Near-Infrared Detector (GROND), a 2.2 m European 
telescope, is capable of observing GRB afterglows simul- 
taneously in 7 wavebands from 400 nm to 2.2 /im, and 
is ideally suited for the type of study we have presented 
in this paper, provided contemporaneous X-ray measure- 
ments can be coordinated. 

The presence of the 2175 A bump can be inferred re- 
gardless of the normalization of the unabsorbed contin- 
uum or its power-law slope. However, these two pieces 
of information are required to derive the shape of the 
extinction curve and the Ay /Ah values, and therein lies 
the power of the IR data. Besides GROND, near-infrared 
(NIR) , ground-based observatories will play a key role in 
such endeavors. 
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Note: Ah is shorthand for A X 10 
O: Light curve fit to May 27.264. 
f: Indicates an upper limit. 
J: Scaled to t = iiR . 



TABLE 1 

Optical/ultraviolet/X-ray Fits & Infrared Observations of GRB 050525A 
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Derived dust extinction in magnitudes 


Band 


Ax (Model A) 
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TABLE 3 

Derived visual extinction per column density 



Model 




A v /Nn (Model A) 


Xr (Model A) 


h 


A v /Nn (Model B) 


Xr (Model B) 


U 






(10^ 22 cm 2 ) 


(N AoC = 6) 




(10 -22 cm 2 ) 


(JVdof = 6) 




MW (R v = 


3.1) 


(4.44 ± 1.41)_ 22 


1.4 


0.8 


(3.36 ± 1.37)_ 22 


1.5 


0.6 


MW {R v = 


4.0) 


(4.10±1.34)_ 22 


1.9 


0.7 


(3.12 ± 1.30)_ 22 


2.1 


0.5 


MW (R v = 


5.5) 


(3.81±1.30)_22 


2.5 


0.5 


(2.92 ± 1.26)_ 22 


2.7 


0.1 


SMC 




(8.55±2.68)_23 


1.1 


1.1 


(6.37±2.52)_ 23 


1.0 


1.1 


LMC 




(1.30±0.40)_22 


0.9 


1.1 


(9.76±3.88)_ 23 


1.0 


0.8 



Note: A h is shorthand for A X 10 b . 
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Fig. 1. — Combined optical and ultraviolet light curves, where data points from the other 6 bands are scaled to the R band model fit 
(dotted curve). Note that the scaling is used only for visualization and was not employed in the fitting algorithm. 
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Fig. 2. — Derived spectral energy distribution (SED) and extinction curve for Model A (see text). Top: SED of the afterglow of GRB 
050525A at t = tin rj 2.3 days after the burst as a function of the observed frequency, v. A power law with an index of /3 = —0.87 ± 0.02 
is fitted to the infrared and X-ray data points. Bottom: Extinction curve, A\, for the host galaxy of GRB 050525A as a function of the 
rest frame frequency, vq. Shown are fits to known dust extinction templates: the Small and Large Magellanic Clouds (SMC and LMC), 
the Milky Way (MW) for different Ry values and a power law fit (with an index of 5 = 1.01 ± 0.31). 
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Fig. 3. — Same as Fig. [2] but for Model B (see text). 



